Ab initio many-body calculations on infinite carbon and boron-nitrogen chains 
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In this paper we report first-principles calculations on 
the ground-state electronic structure of two infinite one- 
dimensional systems: (a) a chain of carbon atoms and (b) 
a chain of alternating boron and nitrogen atoms. Meanfield 
results were obtained using the restricted Hartree-Fock ap- 
proach, while the many-body effects were taken into account 
by second-order M0ller-Plesset perturbation theory and the 
coupled-cluster approach. The calculations were performed 
using 6-31G** basis sets, including the d-type polarization 
functions. Both at the Hartree-Fock (HF) and the correlated 
levels we find that the infinite carbon chain exhibits bond al- 
ternation with alternating single and triple bonds, while the 
boron-nitrogen chain exhibits equidistant bonds. In addition, 
we also performed density-functional-theory-based local den- 
sity approximation (LDA) calculations on the infinite carbon 
chain using the same basis set. Our LDA results, in con- 
tradiction to our HF and correlated results, predict a very 
small bond alternation. Based upon our LDA results for the 
carbon chain, which are in agreement with an earlier LDA 
calculation calculation [ E.J. Bylaska, J.H. Weare, and R. 
Kawai, Phys. Rev. B 58, R7488 (1998).], we conclude that 
the LDA significantly underestimates Peierls distortion. This 
emphasizes that the inclusion of many-particle effects is very 
important for the correct description of Peierls distortion in 
one- dimensional systems. 



I. INTRODUCTION 

Carbon occurs in several allotropic forms depending 
upon the nature of its interatomic bonding. In dia- 
mond it exhibits sp^ hybridization resulting in a three- 
dimensional (3D) structure. In its graphite form it ex- 
hibits sp^ hybridization, resulting in a layered structure. 
While in the carbyne form, the carbon atoms exhibit 
sp hybridization consistent with a predominantly one- 
dimensional (ID), chain-like, character Jd Although, truly 
infinite isolated single chains of carbon atoms do not ex- 
ist in nature, however, there are several experimental ex- 
amples of the tendency of carbon atoms to form ID or 
quasi-lD structures. Lagow et al.0 demonstrated that 
linear chains containing up to 28 atoms can be stabilized 
by adding nonreactive terminal groups. There is strong 



experimental evidence to suggest that carbon chains and 
rings are also precursors to the formation of more cfflftpli- 
cated structures such as fuUerenes and nanotubes.l3"la In 
addition, atom resolved scanning tunneling microscopy 
experiments on the (100) surface oL/3-SiC indicate the 
formation of carbon atomic chainsJj Similarly, simula- 
tions performed on tetrahedral amorphous carbon suig 
faces also indicate the formation of rings and chains .0 
Additionally, for the carbon cluster anions C~ there is 
sufficient experimental evidence to indicate the existence, 
of chains for n < 10, and rings for the larger structures. a 
From a theoretical point of view as well carbon, chains 
and rings have elicited a large amount of interest .ErEJ Let 
us assume an infinitely long carbon chain (C chain here- 
after), with the conjugation direction along the x-axis. If 
the chain were to have equal bond distance (cumulenic 
character with exclusively double bonds), then it will ex- 
hibit metallic character with the occupied band structure 
consisting of: (a) a core band composed of Is orbitals (b) 
a valence a band consisting mainly of sp-hybrids and (c) 
two half- filled degenerate Trpb ands . at , the Fermi level com- 
posed of TTy and TTz 0rbitalsJi3'EBE3O Therefore, the ques- 
tion arises whether such a chain will undergo a Peierls- 
type distortion to become an insulator with bond alterna- 
tion along the chain (acetylenic character with alternat- 
ing single and triple bonds). For a finite ring of carbon 
atoms the same question reads: will it exhibit aromatic 
(equa l hoij ids) or antiaromatic (bond alternated) behav- 
ior ?ml2l'L3 Carbon rings containing AN atoms {N is an in- 
teger) undergo bond alteiaiaikai due to first-order Jahn- 
Teller (JT) distortion.Bililll3'c3 However, the answer for 
a ring containing AN -\- 2 atoms is not straightforward 
because, if it undergoes bond alternation, it can only be 
due to a second order JT effect .lTHIcj Perhaps it is be- 
cause of this uncertainty as to whether the infinite linear 
C chain ring, with large N) will undergo Peierls 

distortion (second-order JT effect) which has ie(i4o nu- 
merous theoretical studies of these systems. □cjEII First 
principles calculations of the infinite C chaip|Jall in two 
categories, namely, the Hartree-Fock (HF)EJ~ri, andjha 
density functional theory (DFT) calculations .BBE3ll3'E£l 
Ab initio HF calculations invariably predict an infinite 
C chain with bond alternation to be energetically pjpf^ 
stable as compared to the equidistant configuration.Eil^o 
However, HF calculations generally have a tendency to 
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overestimate bond alternation which was alsfl-evident in 
our previous studies on tran5-polyacetylene.E3E3 It is only 
the subsequent treatment of electron correlation effects 
which brings the bond-alternation in such systems closer 
to the experimental realityEj Therefore, one would ex- 
pect that the DFT-based approaches, which do include 
electron-correlation effects in an approximate manner, 
would provide an unambiguous picture. However, this 
is not the case. Whereas the local density apprpyiaia^ 
tion (LDA) results of Springborg and collaboratoraljii3'L£l 
place the bond alternation in the range O.lA < S < 0.2A 
(S = Tsingie - rtripie), the reccut LDA study of Bylaska et 
al.M predicts only a very small value i5 < 0.04A for the 
infinite chain. Keeping in view these conflicting results 
obtained for the bond alternation from the DFT calcula- 
tions, and the fact that HF results tend to overestimate 
(5, we believe that a systematic study of the influence 
of electron correlation effects on the bond alternation in 
the infinite C chain is in order. Therefore, we decided 
to apply our ab initio wave-function-based many-body 
methodology uspd. eadier to study a number of conju- 
gated polymers ,Ea'E3L3 to the problem of the infinite C 
chain as well. In addition, to put our correlated results in 
proper perspective, we also performed LDA calculations 
on the C chain using the same basis set which was used 
for the HF and the subsequent many-body calculations. 
These LDA calculations were performed using the CRYS- 
TAL98 programs. Indeed our many-body calculations 
prftdid_A.^« 0.15A close to the findings of Springborg et 
al. Jlj'EaEj bat in complete disagreement with those of By- 
laska et al.M However, ourJJDA calculations performed 
with the CRYSTAL98 codeEj predict a very small bond 
alternation of (5 « 0.04A, in excellent agreement with-the 
LDA value of (5 « 0.0361 reported by Bylaska et al.B 

Another useful concept to quantify the tendency to- 
wards bond alternation is the condensation energy de- 
fined as the difference in energy/atom of the equidis- 
tant configuration and the bond alternating configura- 
tion. Our best value of the condensation energy, ob- 
tained at the many-body level using the coupled-cluster 
approach, is 18.4 milliHartrees/atom. This value is much 
larger than than 0.13 milliHartrees/atom predicted by 
our LDA calculations and also the LDA valuG-©f 0.09 mil- 
liHartrees/atom, reported by Bylaska et alO Thus our 
many-body calculations indicate a much stronger ten- 
dency for the infinite carbon chain to undergo Peierls 
distortion as compared to the LDA calculations. 

In addition to the infinite C chain, in this paper we 
also study an infinite chain composed of boron and ni- 
trogen atoms (BN chain, hereafter). If we consider a 
two atom unit cell, the carbon and the BN chains are 
isoelectronic. Thus if the two chains exhibit different be- 
havior, it will be entirely due to the differences in the 
nature of bonding caused by different atomic species. 
As a matter of fact, starting from well-known carbon- 
based materials, the prospect of obtaining novel mate- 
rial by performing alternate boron and nitrogen sub- 



stitution, has fascinated experimentalists and theoreti- 
cians alike .E3oi23'El Of course, cubic and hexagonal BN, 
which are the counterparts of diamond and graphite, re- 
spectively!,— are well-known materials. In our own ear- 
lier work,E3_as also in the works of Armstrong and col- 
laborators,!^ the BN analogs of trans-polyacetylene and 
polyethylene were investigated theoretically. BN analogs 
of fuUerene-like cages and nanotpbes J,t^ve been studied 
theoretically by several authors.E^O'E^I Experimentally, 
the-J3N nanotubes were first synthesized by Chopra et 
al.j2a while recently, multiwall BN nanotubes, as well as 
fullerene-like BN cages, have been synthesized by Lee et 
al.Jlj in a catalyst-free manner, using rtije laser- ablation 
technique. Additionally, Cote et al.cil calculated the 
optical properties of the BN analogs of photolumines- 
cent conjugated polymers poly(para-phenylene) (PPP) 
and poly(para-phenylenevinylene) (PPV). However, to 
the best of our knowledge, the infinite isolated BN chain 
has not been studied theoretically so far. Therefore, 
in light of the current interest in novel BN compounds 
structurally similar to carbon-based compounds, a par- 
allel study of the infinite BN chain is desirable. Indeed, 
we find that the infinite BN chain has properties dif- 
ferent than the infinite C chain in that it does not ex- 
hibit bond alternation. This result is similar to our ear- 
lier result where we found that the BN analog oLtrans- 
polyacetylene does not exhibit bond alternation.E3 

The remainder of this paper is organized as follows. 
In section ^ the applied methods and computational de- 
tails are briefly described. The results are then presented 
and discussed in section HI. Finally, our conclusions are 



presented in section tV 



II. APPLIED METHODS AND 
COMPUTATIONAL DETAILS. 

For the linear infinite C chain we chose a unit cell con- 
sisting of two carbon atoms, while for the BN chain a 
corresponding unit cell consisting of one boron atom and 
one nitrogen atom was considered. For both systems HF 
calculations were performed first to obtain the results at 
the mean-field level. For the dimerized C chain, as well 
as for the BN chain, the HF calculations were performed 
in the infinite system limit using both the Bloch-jOrbiial- 
based electronic structure program CRYSTAL J^'CJ as 
well as_our own Wannier-function-based program WAN- 
NIER.E3 Both these programs use the linear-combination 
of atomic orbital (LCAO) approach based on Gaussian- 
type basis functions. However, for the equidistant C 
chain, we encountered convergence difficulties with the 
CRYSTAL program when diffuse basis functions were 
used. Therefore, for the metallic C chain, we obtained 
HF results by performing finite cluster calculations using 
clusters of increasing size, and by ensuring that the con- 
vergence with respect to the cluster size was achieved. To 
further check the accuracy of the finite-cluster approach, 
the finite-cluster HF calculations were also repeated for 
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the case of insulating chains for which the infinite system 
HF results (CRYSTAL and WANNIER) were available, 
and excellent agreement between the two sets of results 
was obtained. 

As mentioned earlier, for the purpose of comparison 
with other approaches, we also performed LDA calcula- 
tions on the infinite C chain to study its bond-alternation 
properties. These LDA calcuiations were performed 
with the CRYSTAL98 prograrrS, employing the same 
6-3 IG** basis set as was used in the HF and the sub- 
sequent many-body calculations. Both in the HF and 
the LDA calculations perfcii-mcd with the Bloch-orbital- 
based CRYSTAL program,L§E3 particular attention was 
paid to the convergence of total energies with respect to 
the number of k-points used in the Brillouin-zone (BZ) 
integration. The number of k-points was increased un- 
til the total energy/cell had converged at least up to 
1.0 X 10~^ Hartrees. We noticed that using thirty k- 
points in the irreducible part of the BZ was sufficient 
to achieve aforesaid convergence both for HF and the 
LDA . calc iulations performed with the CRYSTAL pro- 
gram.E3c3 Moreover, the excellent agreement observed 
among the HF energies/cell computed using tiue CRYS- 
TAL program, our own WANNIER programj23 and the 
finite-cluster-based HF calculations, leaves no doubt that 
the convergence in energies has been achieved. 

The many-body calculations beyond HF were per- 
formed only within the finite-cluster model. However, 
since the correlation effects are highly localized in real 
space, convergence with respect to the cluster size is 
rarely a problem here. Moreo"j££x, in our earlier calcula- 
tions on the infinite LiH chainj^j and the bulk LiHjfj we 
had carefully compared the correlated calculations per- 
formed simultaneously on the infinite system and its fi- 
nite cluster, and found excellent agreement between the 
two sets of results. The issue of the convergence of the 
present set of results with respect to the cluster size is 
very important, and, therefore, will be discussed again in 



bonds, was terminated with two H atoms on each end, 
and the energy/cell was computed using the expression 
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section 

In order to perform the correlated calculations, 
three different many-body approaches, viz., second-order 
M0ller-Plesset perturbation theory (MP2), coupled- 
cluster singles and doubles (CCSD), and coupled-cluster 
singles and doubles with perturbative treatment of the 
triples (CCSD(T)) were used. All the finite-cluster-based 
calculations, both at the HF and the correlated level, 
were performed with the MOLPRO molecular orbital 
ab initio program packageo In order to minimize the 
boundary effects associated with the finite clusters, we 
extracted the energy per unit cell from the finite-cluster 
calculations by computing the total energy differences be- 
tween the clusters containing n, and n+1 unit cells, with 
as large a value of n as feasible. In view of the typical ac- 
curacy of the periodic HF codes, i.e., 1 milliHartree per 
atom, we consider the results to be converged with re- 
spect to n, when the (total/correlation) energy per unit 
cell is stable within this accuracy. The undimerized C 
chain, which consists of C atoms connected by double 



E = lim AEn = lim 
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The insulating C chain consisting of alternating single 
and triplet bonds, was terminated with a single CH bond, 
leading to the formula 



E = lim AEn = lim 



E{C2n+2H2) — E{C2nH2) 



(2) 



As far as the BN chain is concerned, its bonding pat- 
tern is not obvious. Therefore, several end-termination 
schemes to saturate the dangling bonds were explored. 
These included terminating the chains with: (a) one H 
atom on each end (b) two H atoms on each end, and (c) 
unequal number of H atoms on the two ends. However, 
converged energies per cell were found to be quite insen- 
sitive to the end geometry of the chain. Therefore, in 
close analogy with the undimerized C chain, we finally 
settled with terminating the BN chain with two H atoms 
on each end, and the expression 



E = lim AEn = lim 

n — ^oo n — >oo 



E{{BN)n+iHi) - E{{BN)nHi) 

(3) 



was used to compute its energy per unit cell. As dis- 
cussed in the next section, the HF energy/cell of the 
BN chain obtained by this method is in excellent agree- 
ment with the ones obtained from the infinite-chain cal- 
culations performed with the CRYSTAL and WANNIER 
programs. Therefore, we believe that the aforesaid ter- 
mination scheme, which effectively models a BN chain as 
composed of double bonds, is a sound one. 



III. RESULTS AND DISCUSSIONS 

Calculations on both the systems were performed with 
6-310** basis sets. For the HF calculations on the infi- 
nite insulating chains, we had to increase the outermost 
p-type exponents of the original basis sets because of the 
linear-dependence related problems. This led to 0.29871 
for carbon (original 0.16871), 0.16675 for boron (origi- 
nal 0.12675), and 0.29203 for nitrogen (original 0.21203). 
The outermost s-type exponents of the original basis set 
were left unaltered. The d-type exponents used for each 
atom were: 0.8 for carbon, 0.6 for boron, and 0.8 also 
for nitrogen. These d exponents were preferred over the 
originaLd-type exponents available from the MOLPRO 
libraryJ13 because they led to lower total energies. For 
the correlated calculations, we used the original basis set 
to compute the energy/cell and the correlation energy 
was computed by subtracting from it the HF energy/cell 
obtained by finite-cluster calculations done also with the 
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original basis set. In this manner, we ensured that no 
spurious contributions to the correlation energy were ob- 
tained because of the use of two different basis sets during 
these calculations. Next we discuss the results obtained 
for the two systems in detail. 



A. Carbon Chain 

Our final results for the equidistant and dimerized C 
chains, employing the wave-function-based approaches 
HF, MP2, CCSD, CCSD(T)), are presented in tables 
_ and H, respectively. Results of our WpA calculations 
performed with CRYSTAL98 programj2j along with the 
comparison of our results with those of other authors, 
are presented in table III . As mentioned earlier we expe- 
rienced convergence difhculties with the CRYSTAL pro- 
gram for the undimerized chain, so HF results for that 
system were obtained by the finite-cluster approach only. 
Moreover, all the correlated calculations were also done 
in a finite-cluster model. Therefore, it is important for us 
to demonstrate the convergence of all the finitc-clustcr- 
based results for both the undimerized, and the dimerized 
C chain, as a function of the cluster size. 

Total HF energy/cell as a function of the number of 
unit cells n in the cluster is plotted in Fig. ||, and it is 
clear that the convergence with respect to the cluster size 
has been achieved both for the dimerized and undimer- 
ized chains, at least, to the milliHartree level. Table|| also 
presents the total energy /cell for the undimerized chain 
including contributions of the electron correlation effects, 
while the corresponding data for the dimerized chain is 
presented in table ||. The convergence of the correlation 
energy/cell computed by the CCSD(T) method for both 
types of chains, with respect to the cluster size, is pre- 
sented in Fig. H. It is clear from these figures that the en- 
ergies for the dimerized chain converge very rapidly, con- 
sistent with our intuitive picture of localized electrons in 
such systems. The convergence in the undimerized chain, 
although comparatively slower, has also been achieved, at 
least, to the milliHartree level. 

The optimized geometries reported in this work were 
obtained by the usual procedure of first performing sev- 
eral total energy calculations for various geometry pa- 
rameters, and then fitting the results by a least-squares 
procedure to polynomials of suitable degrees. The fact 
that the finite-cluster results at the HF level agree fully 
with those obtained by the CRYSTAL and WANNIER 
programs gives us added confidence in their correctness. 
It is clear that the dimerized structure is energetically 
much more stable as compared to the metallic one, at 
ah levels of calculations (SCF, MP2, CCSD, CCSD(T)). 
At the HF level we obtained a percentage bond alter- 
nation <5(%)=14.5% (where S{%) = (S/rav) x 100, with 
rav = irsingie + r triple) / , and at the CCSD(T) level 
of correlation it reduces to 11.8%. The condensation 
energy, defined earlier as the difference in the total en- 



ergy/atom between the optimized undimerized and the 
dimerized C chains, at various levels of correlation treat- 
ments, is obtained to be £^cond(HF) = 8.8 milliHartrees, 
-Eco„d(MP2) = 19.25 milliHartrees, Econd{GCSV)) = 
20.35 milliHartrees, and Econd{GCSV){T)) = 18.4 milli- 
Hartrees. Thus, based on the CCSD(T) approach, which 
is the most sophisticated correlation approach that we 
have used in these calculations, our final result for the 
condensation energy of an infinite C chain of 18.4 mil- 
liHartrees/atom. As reported in table |l|, the value of 
condensation energy obtained from ourXpA calculations 
performed with CRYSTAL98 programN was 0.13 milli- 
Hartrccs/atom, which is very close to the LDA vaiue of 
0.09 milliHartrees/atom reported by Bylaska et al.JlJ and 
much smaller than our many-body results based upon the 
CCSD(T) method. 

The summary of our many-body and LDA results, 
along with those obtained by other authors is presented 
in Tabic [II. It is clear from the table that the results 



obtained by various authors on the optimized geometry 
parameters of the C chain vary quite a bit, however, al- 
most all of them report the bond alternation to be more 
than 10%. The exceptions to this result are (a) our own 
LDA calculatipa (<^(%) = 1-6), (b) LDA calculations of 
Eriksson et al.Ei4|5(%) = 6.0), and (c) LDA calculations 
of Bylaska et al.B((5(%) = 2.9). Thus these recent LDA 
results, which predict rather small bond alternation for 
the C chaiiu-are in contradiction with the older LDA cal- 
culationsli§E3, as well with our many-body results. We 
believe that one should expect a rather small bond al- 
ternation from the LDA based calculations, because the 
homogeneous-electron-gas approximation therein leads 
preferentially to a metallic ground state. At this aoint, 
we would like to mention that recently Torelli et aljj per- 
formed calculations on the optimized ground state ge- 
ometries of C4Ar+2 rings using a quantum Monte Carlo 
method and their results for the bond alternation for 
the Ci8 and C22 rings were 7%. Although, our final 
CCSD(T) result of 11.8% bond ahernation for the in- 
finite C chain is more than that obtained by Torelli et 
al.EI for the large C rings, yet we believe that qualitative 
trends indicated by most of the calculations suggest a sig- 
nificant bond alternation in the infinite C chains. In this 
context we would like to mention that for the finite ring 
systems containing 4A^-|-2 carbon atoms, the Hiickel the- 
ory predicts that the bond alternation due to the second- 
order Jahn- Teller distortion will increase with increasing 
N ^ eventually coalescing with the vaiue corresponding 
to the Peierls distortion, for iV ^ oo.Ea Discounting the 
basis-set effects, in our opinion, our value of bond alter- 
nation of 11.8% for the infinite C chain, as compared to 
7% reported by ToreUi et al.El for the C22 rings, is indica- 
tive that the aforesaid Hiickel result is valid even when 
the many-body effects have been included. 

Other authors have not reported the value of the con- 
densation energies obtained in their calculations. But, 
based on the trends visible from the calculations on bond 
alternation, we believe that it is safe to assume that the 
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condensation energy for the infinite C chain is in ac- 
cess of 10 milliHartrees/atom. Thus for condensation 
energgf, our many-body results, as well that of the Torelli 
et ahlj for the finite rings, are quite different from our 
LDA results of 0.14 milliHartree/atom_and 0.1 milli- 
Hartrec/atom reported by Bylaska et al.EJ 



B. BN Chain 

As mentioned earlier, in nature there appears to be 
a one-to-one correspondence between several structures 
of carbon and BN allotropes such as diamond vs. cubic 
BN, graphite vs. hexagonal BN, and carbon nanotubes 
vs. BN nanotubes etc. Therefore, it will be natural to 
assume that a BN chain, similar to a carbon chain, will 
also have a strictly linear geometry exhibiting sp-type 
hybridization with, 180° bond angles. However, this be- 
ing the first theoretical study of this system, we decided 
to be careful with this assumption. Therefore, in our 
calculations, along with the linear geometry, we also ex- 
plored the energetics of various zig-zag structures of the 
BN chain. However, we found that in all cases the zig-zag 
structure led to higher energies compared to the linear 
one, leading us to the conclusion that if the BN chains 
exist, similar to the C chains, they will have a linear 
structure. Table IV presents the results of our calcu- 



lations for the BN chain both at the HF and the cor- 
related levels. Since we obtained perfect agreement on 
the HF energy/ceU for the BN chain from WANNIER, 
CRYSTAL, and finite-cluster-based approaches, we are 
again very confident of the correctness of our correlated 
studies performed on finite BN clusters saturated with 
hydrogens on the ends. To further confirm this point, 
in Fig. I we present the variation of the correlation en- 
ergy/cell of BN chain, as a function of the number of 
unit cells in the cluster. It is obvious from the figure 
that satisfactory convergence in correlation energy has 
been achieved. At the CCSD(T) level we obtain a cohe- 
sive energy of 11.72 eV/cell, which is somewhat greater 
than the value of 11.05 eV/cell obtained for the dimer- 
ized C chain. This means that from an energetic point 
of view, the BN chain should at least be as stable as an 
isolated C chain, and, therefore, the probability of find- 
ing linear sp-type hybridized structures composed of BN 
units should be as high as that of the C units. 

The other remarkable aspect of the BN chain, as com- 
pared to the C chain, is that it does not exhibit any 
bond alternation. This result is identical to our earlier 
theoretical study on polyiminoborane, which is the BN 
substituted analog of frans-polyacetylene.E3 In contrast 
to frans-polyacetylene which exhibits a bond alternation 
close to 0.1 A, in case of polyiminoborane we obtained 
optimized geometries with-equidistant BN atoms. lj Of 
course, the Peierls theorerrc3 does not apply to BN chain 
or polyiminoborane, because these systems are insulators 
even without bond alternation.cll However, the absence 



of bond alternation in BN compounds can be understood 
from a different line of reasoning as well. BN systems, 
being heteronuclear in nature, have significant ionic con- 
tributions to bonding as opposed to the purely covalent 
carbon-based materials. Therefore, in BN systems in ad- 
dition to the covalent interactions, electrostatic attrac- 
tion between the B and N sites will also exist because of 
a degree of ionicity. It is this electrostatic interaction due 
to ionicity in the BN systems which, in our opinion, leads 
to equal bond lengths. The optimized BN distance at the 
CCSD(T) level was obtained to be 1.3 A, which is slightly 
more than the optimized average bond distance of 1.28 
A obtained for the infinite C chain. For the linear chain, 
our optimized BN distance of 1.3 A should be compared 
to the in plane BN distance of 1.45 A in the hexagonal 
BN, and the BN distance of 1.57 A observed in the cubic 
BN. The significant difference in these three BN distances 
is consistent with the fact that the BN chain, hexagonal 
BN and the cubic BN originate from very different bond- 
ing schemes based on sp, sp^, and sp3 hybridizations, 
respectively. 



IV. CONCLUSIONS 

In conclusion, we have presented an ab initio study of 
the ground state properties of the infinite carbon and 
boron-nitrogen chains. The influence of electron cor- 
relation effects was included by several approaches in- 
cluding the coupled-cluster method accounting for up 
to triple excitation operators. As far as the carbon 
chain is concerned, all our many-body calculations pre- 
dict a dimerized ground state, consistent with the phe- 
nomenon of Peierls distortion. Thus our calculations 
on the inflnite C chain based upon many-body meth- 
ods such as the coupled-cluster approach arc in agre&j 
ment with the Monte Carlo calculations of Torelli et al.EI 
performed on carbon rings, which also revealed a strong 
tendency of these systems to exhibit bond alternation 
in the infinite chain limit. We also performed a parallel 
LDA calculation on the infinite C chain which predicted 
very weak dimerization, in excellent quantitative agree^ 
ment with the recent LDA results of Bylaska et al.EJ 
This, in our opinion, establishes in most unambiguous 
manner that LDA underestimates Peierls distortion in 
one-dimensional systems, and that powerful many-body 
methods such as the couple-cluster approach are essen- 
tial for the correct description of the broken-symmetry 
ground state. 

In addition, we also presented first principles correlated 
calculations on the ground state of an infinite boron- 
nitrogen chain. Based upon our calculations we conclude 
that the boron-nitrogen chain has a higher cohesive en- 
ergy/cell as compared to the carbon chain, thus making 
BN-based ID-structures energetically at least as favor- 
able as the carbon-based ones. Of course, it is a dif- 
ficult task to synthesize isolated infinite atomic chains. 
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However, similar to the case of carbon, it should cer- 
tainly be possible to test in the laboratory whether fi- 
nite one-dimensional clusters based on BN units exist. 
Along the same lines, it will also be worthwhile to the- 
oretically optimize the geometries of the ground states 
of various BN-based finite clusters. In addition to the 
ground state, theoretical studies of their optically excited 
states will also be useful because they, in conjunction 
with photoabsorption-based experiments, can be used to 
characterize various clusters. We will investigate the elec- 
tronic structure of BN clusters in a future pubhcation. 
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FIG. 1. Total HF energy per unit cell for both the undimer- 
ized (circles) and dimerized (squares) carbon chains computed 
by the finite-cluster approach (cf. Eqs. (|^), and (^), plotted 
as a function of the number of unit cells n. Two carbon atoms 
were assumed per unit cell. The bond distances in these cal- 
culations correspond to the optimized geometries reported in 
tables I and 0. 

A. Abdurahman et al.: Fig 1 
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FIG. 2. Correlation energy per unit cell computed by the 
CCSD(T) approach both for undimerized (circles) and dimer- 
ized (squares) finite carbon chains, plotted as a function of 
the number of unit cells n in the finite cluster. For the 
undimerized chain rdouUe = 1.2QA, and for the dimerized 
chain rsingie = 1.37A, and rtripie = 1.175A. 

A. Abdurahman et al.: Fig 2 
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FIG. 3. Correlation energy per unit cell of finite 
boron-nitrogen chains computed by the CCSD(T) approach, 
plotted as a function of the number of unit cells n in the finite 
cluster. The bond distance was taken to he vbn = 1.301 A. 



A. Abdurahman et al.: Fig 3 
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TABLE I. Total energy Etot (Hartrcc), cohesive energy 
AEcoh (eV) per unit (C — C) and bond lengths(A) for the 
undimerized structure of the infinite carbon chain. Cohesive 
energies were computed with respect to separated atoms. 



Method 




^Eco/i 


^double 


Finite cluster SCF 


-75.5644 


5.733 


1.251 


MP2 


-75.7847 


10.327 


1.268 


CCSD 


-75.7813 


9.570 


1.267 


CCSD(T) 


-75.7996 


10.052 


1.272 
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TABLE II. Total energy Etot (Hartree), cohesive energy 
AEcoh (eV) per unit (C — C) and bond lengths(A) for the 
dimerized structure of the infinite carbon chain. Cohesive 
energies were computed with respect to separated atoms. 



Method 


Etot 


AEcoh 


1^ single 


'^triple 


WANNIER SCF 


-75.5807 


6.176 


1.367 


1.173 


CRYSTAL SCF 


-75.5798 


6.150 


1.369 


1.173 


Finite cluster SCF 


-75.5810 


6.184 


1.360 


1.174 


MP2 


-75.8232 


11.375 


1.337 


1.217 


CCSD 


-75.8220 


10.678 


1.362 


1.197 


CCSD(T) 


-75.8364 


11.053 


1.358 


1.207 



Correlation contributions added to CRYSTAL SCF en- 
ergies. 



TABLE III. Infinite carbon chain: comparison of the 
present results with those of other authors. The per- 
centage bond alternation = (S/rav) x 100, where 

^ ^single ^triple^ aud Tav ij' single Ttriple) l'^- 



Author 



Method TsingU rtriple S{%) Econd 



(A) (A) 



(mHartrees / atom) 



This work 
This work 
This work 
Bylaska et al, ^' 
Eriksson et al 
Springborgrtt al 
Springb 
KarpfenE3 
Kertesz et al.— 
Kertesz et al. 'S, 
Teramae et aBIl, 



RHF 

CCSD(T) 

LDA 

LDA 

LDA 

LDA 

LDA 

RHF 

RHF 

UHF 

RHF 



1.360 

1.358 

1.286 

1.288" 

1.376 

1.503 

1.439 

1.363 

1.145^= 

1.335 

1.339 



1.174 
1.207 
1.246 



14.5 
11.8 
1.6 



1.252" 2.9 



7.80 
18.40 
0.13 
0.09* 



1.296 

1.259 

1.249 

1.198 

1.405^= 

1.185 

1.166 



6.0 

18.0 

14.1 

12.9 

20.4 

11.9 

14.1 



° Calculated from the average bond distance of 1.270 A 
and 6 = 2.9 %. ,-, 
* Upper limit (Ref.B) 

Calculated from the average bond distance of 1.275 A 
and 6 = 20.4 %. 
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TABLE IV. Total energy Etot (Hartree), cohesive energy 
AEcoh (eV) per unit (BN) and B-N bond length tbn (A) 
of the infinite boron-nitrogen chain. Cohesive energies were 
computed with respect to separated atoms. 



Method 


Etot 






WANNIER SCF 


-79.1801 


7.5845 


1.287 


CRYSTAL SCF 


-79.1794 


7.5656 


1.288 


Finite cluster SCF 


-79.1773 


7.5076 


1.287 


MP2 


-79.3968 


11.9478 


1.300 


CCSD 


-79.4017 


11.4723 


1.299 


CCSD(T) 


-79.111:! 


11.7201 


1..301 



Correlation contributions added to CRYSTAL SCF en- 
ergies. 
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